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ABSTRACT 

We present new constraints on the evolution of the early-type galaxy color-magnitude relation 
(CMR) based on deep near-infrared imaging of a galaxy protocluster at z = 2.16 obtained using 
NICMOS on-board the Hubble Space Telescope. This field contains a spectroscopically confirmed 
space-overdensity of Lyman-a and H-a emitting galaxies which surrounds the powerful radio galaxy 
MRC 1138-262. Using these NICMOS data we identify a significant surface-overdensity (= 6.2x) 
of red Jno — i?i60 galaxies in the color-magnitude diagram (when compared with deep NICMOS 
imaging from the HDF-N and UDF). The optical-NIR colors of these prospective red-sequence galaxies 
indicate the presence of on-going dust-obscured star- formation or recently formed (< 1.5 Gyr) stellar 
populations in a majority of the red galaxies. We measure the slope and intrinsic scatter of the CMR 
for three different red galaxy samples selected by a wide color cut, and using photometric redshifts 
both with and without restrictions on rest-frame optical morphology. In all three cases both the rest- 
frame U — B slope and intrinsic color scatter are considerably higher than corresponding values for 
lower redshift galaxy clusters. These results suggest that while some relatively quiescent galaxies do 
exist in this protocluster both the majority of the galaxy population and hence the color-magnitude 
relation are still in the process of forming, as expected. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: high-redshift — galaxies: 
stellar content — galaxies: clusters: individual: MRC 1138-262 



1. INTRODUCTION 

The color-magnitude diagram is a powerful diagnos- 
tic of galaxy evolution and formation. The presence, 
as early as z ~ 1.5, of a prominent and low-scatter, 
'red-sequence' (RS) in galaxy clusters places useful con- 
straints on the possible evolutionary pathways in galaxy 
color and luminosity ( Mullis et al.1 120 05: Sta nford et al.l 
[20051 [20061: iBell et al.ll2004t iFaber et"al"ll2007f) . The red 
colors of the primarily early-type RS galaxies are due to 
the observed filters spanning the 4000A spectral break. 
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The universality and prominence of the RS in appro- 
priately chosen filters have been used to di scover high- 
redshift clusters (e.g.. Idadders fc Yeell2005l ). Moreover, 
the defining characteristic of galaxy clusters, i.e., the 
large numbers of galaxies all at the same redshift, allows 
the slope and intrinsic scatter of the RS to be measured 
with great precision. Based on studies of galaxy clus- 
ters at z < 1.3, the slope of the RS does not appear to 
evolve and therefore is more likely the by-product of the 
mass-met allicity relation as obs erved in local galaxy sam- 
ples (e.g.. iTremonti et al.ll2004T ) rather than the result of 
a mass-age trend. The scatter, however, is likely due 
to the fr actional age difference s between the RS galax- 
ies (e.g.. [Blakeslee et al.ll2003l ). By constructing a set 
of model galaxies with different star-formation histories 
and timescales it is possible to fit for the mean epoch of 
last significant star-formation by matching the intrinsic 
scatter of the RS. Such studies at z ~ 1 h ave derived for- 
matio n redshifts of Zf nr m ~ 2.0 — 2.5 (e.g..lBlakeslee et al.1 
120061: Ivan Dokkum fc van der Marell2007t ). At redshifts 
beyond z ~ 1.5, however, the 4000A-break moves into 
the near-infrared and galaxy clusters, and therefore the 
RS, have not been observed closer to the inferred epoch 
of formation for early-type galaxies. Hence, to uncover 
the younger or forming red-sequence at higher redshifts 
requires deep near-infrared imaging of suspected (or, 
preferably, confirmed) protocluster fields. 

We have undertaken a NICMOS imaging program to 
study the red galaxy population in a protocluster at z = 
2.16. Broad and narrow-band imaging, both in the opti- 
cal and near-infrared, of the field surrounding the power- 
ful radio galaxy MRC 1138-262 (z = 2.16) have identified 
more than 100 candidate companion galaxies. This tar- 
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Fig. 1. — Jno — -ffieo vs. Hieo color-magnitude diagram for the MRC 1138-262 NICMOS field (open black circles). The large yellow star 
is the radio galaxy itself. The blue background points are from the NICMOS data of the Ultra Deep Field and Hubble Deep Field North. 
The deep fields cover 2.5x the area of our observations. Also shown are the spectroscopically confirmed Ha (orange filled circles) and 
Lyman-o (green filled circles) emitters. The three dot-dash lines show where the CMRs of lower redshift clusters would lie under different 
assumptions. The top line is the Coma cluster with no evolution, simply redshifted to z = 2.16 and observed through the NICMOS filters. 
The next line down is the same but for the z = 1.24 cluster RDCS1252. Finally, if we passively de-evolve RDCS1252 to redshift z = 2.16 
(almost exactly two Gyr), assuming a median age for the 1252 galaxies of about 3 Gyr (or Zf orm ~ 3), we get the third line. 



get served as the proof-of-concept for the successful VLT 
Large Program summarized in Venemans et al. (2007). 
There are surface-overdensities of both line-emitting can- 
didates (Lyman-a and Ha), X-ray point sources, sub- mm 
selected galaxies and r ed optical-near-infrared galaxi es 
dPentericci et al 1 120021: iKurkl 120031: iKurk et al.ll2004aH bl: 
iCroft et a l. 2005; S tevens et aI1l2003h . Fifteen of the Lja 
and 9 of the Ha emitters have now been spectroscop- 
ically confirmed to lie at the same redshift as the ra- 
dio galaxy. The / — i^-selected Extremely Red Objects 
(EROs; I — K > 4.3 Vega) seem concentrated around 
the RG but have no spectroscopic redshifts at this time. 
However, by obtaining deep images through the NIC- 
MOS Jno and Hiqo filters, which effectively span the 
4000A-break at z = 2.16, accurate and precise colors 
and basic morphological parameters can be measured for 



the red galaxy population. In this paper we present the 
first results from this project. The article is organized 
as follows: in Section §?? we describe the data and their 
reductions, in Section §?? we present the comparison be- 
tween the red galaxy counts in this field and in deep field 
data, in Section §?? we present the full color-magnitude 
diagram and our fits to the "red sequence." We use a 
(0 A) ^m) = (0.73,0.27), H = 71 km s" 1 Mpc" 1 cos- 
mology throughout. At z = 2.16 one arcsecond is equiv- 
alent to 8.4 kpc All magnitudes are referenced to the 
AB system (|Okelll97"4T ) unless otherwise noted. 

2. OBSERVATIONS, DATA REDUCTIONS AND 
PHOTOMETRY 

The NICMOS instrument on-board HST is capable 
of deep near-infrared imaging more quickly than from 



The Nascent Red- Sequence at z ~ 2 



3 




1 



' J 1 10 n 160 

Fig. 2. — 9435 — ^814 vs - JllO — Hi60 color-color diagram us- 
ing the ACS and NICMOS data. Arrows represent limits where 
the galaxy is only detected in a single band for that color. Filled 
circles indicate spectroscopically-confirmed Lyman-a (green) and 
H-a (orange) emitting protocluster members. The yellow star is the 
radio galaxy. The blue, green and red grids indicate the regions oc- 
cupied by galaxies with an exponentially-decaying star-formation 
rate r = 0.15 Gyr (red), r = 0.4 Gyr (green) and a r = 1000 
Gyr (blue) at z = 2.16 for three ages (0.1, 1 and 3 Gyr) and two 
different extinctions (E(B - V) = 0.0, 1.0). 

the ground but with a relatively small field-of-view 
(51" x 51"). In the case of MRC 1138-262, we know 
that gala xies are overdense on the scale of a few ar- 
cminutes (|Kurk et alj|2004at ICroft et all [20051) and are 
thus well-suited for observations with NICMOS camera 
3 on HST. We used 30 orbits of HST time to image 
10 of the 24 confirmed members and ~ 70 of the can- 
didate (narrow-band excess sources and EROs) proto- 
cluster members in both the J%xo and Hi 60 filters. We 
used seven pointings of NICMOS camera 3 in both fil- 
ters and one additional pointing in Hiqq alone. This 
single 'outrigger' -ffi60 pointing was included to obtain 
rest-frame optical morphological information for a small 
concentration of candidate members. These observations 
reach an AB limiting magnitude (mw a ; 10a, 0"5 diam- 
eter circular aperture) of toio ct = 24.9 mag in Jno and 
Wio<t = 25.1 mag in i/i6o- The same field was imaged in 
the (7475 (toioo- = 27.5 mag) and I^u {mw<r = 26.8 mag) 
filters using the Wide-Field Channel of the Advanced 
Camera for Surveys on HST as part of a Guaranteed 
Time program (# 10327; Miley et al. 2006). These opti- 
cal data are useful for their higher angular resolution and 
their coverage of the rest-frame far-UV, thus extending 
the observed SEDs of candidate protocluster members to 
shorter wavelengths where young stars and on-going star- 
formation dominate the emitted spectrum. In particular, 
the (7475 and /si4 data allow us to partially differentiate 
obscured star-formation from evolved stellar populations 
in the candidate RS galaxies. 

The NICMOS images were reduced using the on-the-fly 
reductions from the HST archive, the IRAF task PED- 
SKY and the dither/drizzle package to combine the im- 
ages in a mosaic. The dither offsets were calculated using 
image cross-correlation and were refined with one further 
iteration of cross-correlation. Alignment of the pointings 



relative to each other was accomplished using a rebinncd 
version of the ACS /si4 image as a reference. The final 
mosaic has a pixel scale of O'.'l. Galaxies were selected 
usin g the .Hipp-band image f or detection within SExtrac- 
tor (|Bertin fc Arnoutsl [19961. We used a 2.2c detection 
threshold with a minimum connected area of 10 pixels. 
We also corrected the NI CMOS data fo r the count-rate 
dependent non-linearity (|de Jond I2006D . Total galaxy 
magnitudes were estimated by using the MAG_AUTO 
values from SExtractor. 

The Jno — J?i6o colors were det ermined by running 
SExtractor (jBertin &; Arnouti [l996) in two-image mode 
using the -£/i60 image for object detection and isophotal 
apertures. The Jno image was PSF-matched to the i/i60 
band. The resulting colors and magnitudes are shown 
in Figure [TJ For galaxies which are not detected at 2a 
significance in the Jno-band (those to the right of the 
thick dashed line, representing Jno, tot > 26.7, in Fig. [I]) 
we consider the color to be a lower limit. 

We also measured similarly PSF-matched, isophotal 
colors using the two ACS bands and have used them to 
construct a (7475 — Isu versus Jno — -Hi60 color-color di- 
agram (Figure [2]). We compared these colors to those 
of model SEDs for different ages, star-formation his- 
tories and dust extinctions. Using the 2007 Chariot 
& Bruzual population synthesis models we have con- 
structed spectral energy distributions for galaxies with 
an exponentially-decaying star-formation rate with time 
constants of r = 0.15,0.4,1000.0 Gyr (the red, green 
and blue grids in Fig. [2] respectively). Each model's 
colors are calculated for ages of 0.1, 1 and 3 Gyr and 
for E(B — V) — 0.0 and 1.0. Aging of the population 
moves primarily the Jno — -Hi60 color to the red while 
the dust extinction significantly reddens the (7475 — Jsi4 
color. From this analysis it appears that the r = 0.4 Gyr 
model represents well the colors of a majority of the red 
Jno - #160 galaxies. 

To extend the wavelength coverage for the protoclus- 
ter galaxies we also incorporated ground-based J7„-band 
data from LRIS-B on the Keck telescope, if s -band imag- 
ing from VLT/ISAAC and three band IRAC imaging 
(the 3.6, 4.5 and 5.8 /im bands) from the Spitzer Space 
Telescope. The Keck £/-band data (PI W. van Breugel) 
were obtained in late January and early February of 2003. 
The ISAAC data (PI G. Miley) were taken in Period 73 
in service mode. The Spitzer data are from the IRAC 
Guaranteed Time program (PI G. Fazio, Program #17). 
We have smoothed the imaging data for all bands, apart 
from the IRAC data, to match the resolution of the U n - 
band image (approximated by a FWHM ~ 1" Gaussian). 
We then used SExtractor to measure galaxy magnitudes 
within a C'5 radius circular aperture for each smoothed 
image. To incorporate the IRAC data, which has much 
poorer angular resolution, we derived aperture magni- 
tudes which were then corrected to match the smoothed 
data. These aperture corrections were derived using the 
photometric curves-of-growth for 20 stars in the field. 
The resulting catalog was used to generate photometric 
redshift estimates as described below in Section § ??. 

3. PHOTOMETRIC REDSHIFTS 

We have used the ACS (3475, hu), NICMOS (J uo , 
f/i6o), ground-based t/„-band from Keck/LRIS-B, K$- 
band imaging from VLT/ISAAC and Spitzer/IRAC 
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Fig. 3. — upper panel: Distribution of high-confidence (> 95%) 
photometric redshifts and their selection function assuming a uni- 
form N(z) for our model template galaxies (yellow curve) for 
1.1 < (Jno - H 160 ) < 2.1 galaxies in the MRC 1138-262 NIC- 
MOS field. The peak between z = 2.1 and z = 2.4 is statistically 
highly significant, middle panel: Sum of redshift probability distri- 
butions for all the galaxies in the upper panel. 38.5% of the total 
probability is contained in the redshift interval from 2 to 2.3. lower 
2 panels: two examples of the probability distribution function for 
individual galaxies 

imaging to estimate photometric redshifts for our H\qo- 
band selected sample. We input a catalog of aperture 
galaxy magnitudes, based on the matched, smoothed 
images described above, into the Bayesian photomet- 
ric redshift code (BPZ) of Bcnftez (2000) using a uni- 
form prior. We felt that the default prior, based on 
optical galaxy selection and spectroscopy in the HDF- 
N, would not necessarily represent the redshift dis- 
tribution for our near-infrared selected galaxies. We 
generated our own extensive set of template spectral 
energy distributions using the models of Chariot & 
Bruzual (2007). All these SEDs are r models with 
values for r = [0.15,0.4,1.0,2.0,1000.0] Gyr and ages 
= [0.05, 0.1, 0.5, 1.0, 2.0, 3.0] Gyr. We also included mod- 
els with internal dust extinction ranging from E(B — 
V) = [0.0,0.1,0.3,0.5,0.75,1.0] mag and metallicity of 
(Z/Zq) = [0.3,1.0,2.5]. We focused particular attention 
on the Jno — -Hi60 selected surface-overdensity. In the 
upper panel of Figure [3] we present the high confidence 
(> 95%) photo- z distribution for the NIR-color selected 
(1.1 < (Jno — Hiso) < 2.1) subsample. We ran extensive 
simulations by redshifting our template set, adding ap- 
propriate photometric errors and using BPZ to recover 
the redshifts. The yellow curve represents the redshift 
selection function for this color cut, template set and fil- 
ters assuming that these model galaxies follow a uniform 
N(z) over this redshift interval. The simulation results 
were free of significant systematic errors and the random 
errors are estimated to be Sz/z ~ 0.1. Based on these 
SED fits, the approximate luminosity-weighted ages of 
the red galaxies lie between 1 and 2.5 Gyrs and their 
stellar masses are of order a few xl0 10 M Q . These stel- 
lar masses are reasonable as are the absolute magnitudes 
(see Figure O- More detailed SED modeling is deferred 
to a future paper. 



There is a clear excess of galaxies between z = 2.0 and 
z = 2.5. For each galaxy fit by BPZ we have generated 
the full redshift probability distribution. In the lower 
panel of Fig. [3] we show the 7Ji6o-band weighted-average 
of these probability distributions. There is a clear peak 
(containing 38.5% of the total probability compared to 
only 17% of the total selection function in the same red- 
shift interval) between z = 2.0 and z = 2.3, consistent 
with the significant peak in the redshift histogram itself. 

4. NICMOS GALAXY MORPHOLOGIES 

NICMOS camera 3 provides good angular resolution 
over its entire field-of-view. The FWHM of the PSF in 
our final mosaic is w 0V27. T o exploit this reso lution we 
have used the GALFIT code (|Peng et al.ll2002fl to fit an- 
alytic Sersic surface-brightness profiles (|Sersidll968T ) to 
all the i?i6o < 24.5 sources in our iJi 60 -band mosaic. A 
model point-spread function was created for each of these 
galax ies individu ally by generating a TinyTim simulated 
PSF (Krist 1993) at the galaxies' positions in each expo- 
sure and then drizzling these PSFs together in exactly the 
same fashion as for the data themselves (see Zirm et al. 
2007). We restricted the Sersic index, n, to be between 1 
and 5. We will present a full analysis of the morphologies 
of these galaxies in a future paper. For the current work, 
we use these derived sizes and profile shapes to assist us 
in selecting the morphological "early-type" members of 
the red galaxy population. 

In Figure |4] we show the distribution of galaxy ages de- 
rived via these SED fits as parametrized by the r value 
for the best-fitting model for those galaxies with high 
and low Sersic index (n > 2.5, red line, and n < 2.5, blue 
line) . It is clear that while there is substantial overlap be- 
tween these distributions they are not identical and that 
they differ in the sense that one might expect, namely, 
that the concentrated galaxies appear to be comprised of 
older stellar populations. This trend gives us some con- 
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Fig. 4. — Distribution of inferred stellar ages (in terms of r) 
for both the concentrated (n > 2.5, red line) and diffuse (n < 2.5, 
blue line) galaxies which are well-resolved in the NICMOS data. 
Constant star-formation models, for which the e-folding time is 
infinite, are placed at the left-hand edge of the plot. The blue 
and red distributions are quite different. Of particular note is that 
the most evolved galaxies generally have high n while the low n 
galaxies dominate the star-forming population. 
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Fig. 5. — Histogram of the color distributions for the 1138 and deep fields (blue). The deep field data has been normalized by total area 
to the 1138 data. Note the clear excess of red galaxies in the 1138 field. At 1.1 < (Jno — J^i6o) < 2.1 (horizontal dotted lines) for galaxies 
brighter than the 2a Jno-band limit (dashed line) there is an overdensity of a factor of 6.2 in the 1138 field. 



fidence in trying to select the "early-type" galaxies using 
these data which is important for our discussion of the 
color-magnitude relation in Section ??. 

5. SURFACE OVERDENSITY OF RED GALAXIES 

To compare this protocluster field to more generic 
'blank' field data we have compiled catalogs for the 
public NICMOS data in both the Hubble Deep Field 
North (HDF-N) and the Ultra Deep Field (UDF). Fig- 
ure [T] shows the Jno — -Hi60 color-magnitude diagram 
(open black circles) and the color distributions for both 
the MRC 1138-262 and the combined HDF-N and UDF 
galaxy catalogs (blue circles). The deep field data were 
also £/i6o-band selected. The area of the two deep fields 
is roughly 2.5 times the area of our protocluster ob- 
servations. We have applied no correction to the deep 
field number counts to account for clustering in those 
fields. The color histogram in Figure \5\ shows the area- 
normalized galaxy counts from the two deep fields (blue 



line) and from the 1138 field to the same (2cr) limiting 
magnitude of Jno = 26.7 mag (AB). The red dashed line 
shows the difference between the two color distributions. 
It is clear that the radio galaxy field is overdense in red 
galaxies by a large factor. For sources with colors be- 
tween 1.1 < (Jno — Hteo) < 2.1, the horizontal(vertical) 
dotted lines in FigtTUSJ) , and brighter than our Jn -band 
2cr limit (26.7), we calculate an area-normalized overden- 
sity of 6.2 when compared to the deep fields data. We 
note that the exact value of the measured overdensity is 
rather sensitive to systematic color offsets between the 
protocluster and deep field data. A redward shift of 0.05 
for the deep field galaxies would lower the measured over- 
density to 5.0. However, we are confident that these sys- 
tematic offsets remain small (< 0.05 mag) since we have 
used the same instrument, filters, selection technique and 
photometric code with very similar input parameters for 
both the deep field and 1138 datasets. Looking back at 
Figure [2] we can see that many of the spectroscopically- 
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Fig. 6. — Linear fits to the rest-frame U — B (Vega) color-magnitude diagrams for three different sub-samples of the H\ao-hai\d 
selected N1CMOS sample. Panel (A) shows the fit (solid line) and intrinsic scatter (lcr; dotted lines) for a sample selected to have 
1.1 < (Jno — Hi6o) < 2.1. The crossed-out point are those which are rejected as outliers in more than half of the realizations (see §??). 
Both the observed and intrinsic scatter are smaller than the initial color cut. Panel (B) shows the fit and intrinsic scatter for a photometric 
redshift selected sample with 2.0 < z p hot < 2.5. The stars indicate galaxies whose preferred photometric template has an age < 4t, while 
circles represent galaxies with older than this limit. Panel (C) shows the fit and scatter for those galaxies which meet the same redshift cut 
but also are well-resolved with a high Sersic index (n > 2.5) and best-fit by an age > 4r template. 



confirmed line emitters (filled blue circles) and red galax- 
ies in the overdensity are well-represented by the r = 0.4 
Gyr model (green lines) at different ages and extinctions. 

This current work is not the first to observe red galax- 
ies in this field. Kurk et al. (2004) identified a small 
(~ 1.5x) surface overdensity of extremely red objects 
(EROs; I — K > 4.3 Vega magnitudes) using ground- 
based I and K band data. Many of these EROs are also 
identified as red in the NICMOS Jno — #160 color. More 
recently Kodama et al. (2007) observed this field us- 
ing the wide-field NIR imager, MOIRCS, on the Subaru 
telescope. These authors found several bright (presum- 
ably massive) red galaxies over a wider field-of-view but 
to shallower depths than the NICMOS data presented 
here. 24 of their color-selected protocluster candidates 
are within our NICMOS mosaic. 23 of the 24 are iden- 
tified in our data as being red in Jno — iii60- Further- 
more, 18 of the 94 galaxies which satisfy our color criteria 
(and have Jno < 26.67) are also identified by Kodama 
et al. as protocluster candidates. The much larger num- 
ber of red galaxies in the NICMOS data is primarily due 
to fainter galaxies detected at high significance in our 
deeper data. 

6. AN EMERGENT RED-SEQUENCE? 

To study the colors and magnitudes of these galaxies 
in more detail and to possibly identify a red-sequence 
in the 1138 field we have split the galaxies into three 
sub-samples defined by Jno — -Hi60 color, photometric 
redshift and morphology (Sersic index). The first sam- 
ple (Sample A) comprises all 56 galaxies with 1.1 < 
( Jno — -Hi6o) < 2.1 and H wo < 24.5 and includes the ra- 
dio galaxy itself. Sample B is made up of all 28 galaxies 
with a robust photometric redshift between 2.0 and 2.5 
and Jno - H im > 0.75 and H wa < 26.0. This liberal 
color cut is included to select galaxies which comprise 
the large observed surface-overdensity. Finally, sample C 
contains seven galaxies with the same photometric red- 
shift cut but which also have well-resolved f/i60-band 
surface-brightness profiles with Sersic index n > 2.5. All 
of these galaxies' SEDs are also best-fit by models with 
relatively little on-going star-formation. We use a limit 



of (age > 4 x t, cf. Grazian et al 2007 ). Therefore, 
sample C mimics the color, morphological and photo- 
metric redshift selection of early-type galaxies in clusters 
at z < 1. The photometry, photo-zs and sizes of the 
sample C galaxies are listed in Table [TJ their rest- frame 
color-magnitude distributions are shown in Figure [6] and 
the two-dimensional spatial distribution of the Sample 
A galaxies is plotted in Figure [7] We note that because 
the measured overdensity is a factor of 6, we statistically 
expect one of every seven sample A galaxies to be a field 
galaxy. However, this should not effect our results sig- 
nificantly. 

For these three sample selections we have fit a line and 
measured the intrinsic scatter about that best-fit line (see 
Fig. [5]). For comparison to lower redshift galaxy clusters 
we have transformed our observed Jno — -Hi60 color and 
-ffi6o magnitudes into rest-frame U — B and B (Vega), 
respectively, using the following expressions: 

(U - B) Icst = 0.539 x (J 110 - H 16() ) ohs - 0.653 (1) 



M BtTest = i/i 60 ,obs -0. 170 x( J 110 -^i6o)obs -43.625 (2) 

The small color corrections used in these relations were 
derived using a family of r-models with a range of ages 
(0.1-12 Gyr), r (0.1-5 Gyr) and three metallicities (0.4, 1 
and 2.5Z ). 

To fit the "CMR" we used a bootstrap re-sampling 
technique to estimate the error on the fitted slope. Then, 
by assuming that all the red galaxies lie on this fit line, we 
ran Monte Carlo realizations of the contribution of the 
photometric errors to the observed color scatter about 
the fit line, i.e., by fixing a color-magnitude relation 
we calculate the measurement scatter with zero intrin- 
sic scatter. We then calculate the intrinsic scatter by 
subtracting (in quadrature) the estimated measurement 
scatter from the observed scatter. We show these fits 
(solid line) and the intrinsic scatters (dotted lines) for 
the three samples (A, B and C) in Figure [5] The fits 
to both Sample's A and B have nearly identical rest- 
frame U — B slopes, 0.027 and 0.026 respectively, and 
intrinsic scatters (0.10 and 0.12). While these slopes 
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Fig. 7. — Spatial distribution of galaxies (red circles) relative to the radio galaxy MRC 1138-262 (yellow star) for the color-selected 
sample defined in the text (A) and shown in the first panel of Fig. \E\ The irregular black outline encloses the coverage of the 7 NICMOS 
camera 3 pointing with both Jno and ffi60 imaging. 



are comparable to those found for the well-populated 
lower redshift cluster CMRs, the intrinsic scatters are 
considerably higher. However, the scatter measured for 
the eight galaxy Sample C is comparable to that of the 
lower redshift samples but with a much steeper slope 
(0.130). When these scatters are compared to model 
predictions based on lower redshift clusters (specifically 
RDCS 1252.9-2927 at z = 1.24; Gobat et al. 2008) we 
find that the 1138 protocluster has lower than predicted 
scatter. This may suggest that the 1138 protocluster is 
in a more advanced evolutionary state than RDCS 1252 
was at z = 2.2. 

We have calculated three representative color- 
magnitude relations for comparison to the colors and 
magnitudes of the red galaxies (three dot-dash lines in 
Fig- El)- We have taken two lower redshift clusters, Coma 
at z = 0.023 and RDCS 1252.9-2927 at z = 1.24, and 
transformed them to the observed filters and z = 2.16 
under the assumption that the colors do not evolve. In 
this no evolution case (the two dot-dash lines in Fig. [l} 
the CMRs appear at the red edge of the observed over- 
density. There is almost exactly 2 Gyr of cosmic time 
between z = 2.16 and z = 1 .24 in our adopted cos- 
mology. From Blakesle e et al.l ([2003) we know that the 
median redshift of last significant star-formation for the 
RDCS 1252 galaxies is between z = 2.7 - 3.6. There- 
fore, if we observe those galaxies at z = 2.16 they will 
be significantly younger and hence bluer. In fact, this 
passively de-evolved line (bluest dot-dash line in Fig. [TJ 
labeled 'zf orm ~ 3') does fall blueward of the red galaxy 
overdensity. We discuss the implications for these com- 
parisons in Section ??. 

We have also translated the Kodama et al. ground- 
based J — K colors to our NICMOS filters assuming all 
the red galaxies lie at z = 2.16. These bright galax- 
ies also fall along the passively de-evolved line with 
the radio galaxy. We have used our suite of SED 
models to estimate the color transformation from their 
ground-based J — K s to our NICMOS J\ W — Hi m color. 
Roughly, the Kodama et al. bright red galaxies fall 



where the RDCS 1252 passive line crosses our color cut 
at Jno — i?i60 = 1-1. This result hints at a possible 
bi-modality in the red galaxy population of this proto- 
cluster. Namely, that there are faint red galaxies that 
are inconsistent with passively-evolving cluster members 
either due to large amounts of dust, or due to higher 
rcdshifts of formation but that the more luminous proto- 
cluster members may have already finished forming and 
seem consistent with passive evolution to the present-day. 

7. DISCUSSION 

We have identified a (6.2 x) surface-overdensity and a 
corresponding photometric redshift 'spike' of red Jno — 
Hiqq galaxies which are likely associated with a known 
protocluster at z — 2.16. The optical-NIR spectral 
energy distributions of these sources suggest that they 
comprise both evolved galaxies as well as dust-obscured 
star-forming galaxies. Based on our SED fits from the 
photometric redshift determinations, the approximate 
luminosity-weighted ages of these sources lie between 1 
and 2.5 Gyrs and their stellar masses are of order a few 
xl0 10 M Q . Detailed modeling of the SEDs for the pro- 
tocluster population, along with their morphologies, is 
reserved for a future paper. 

Comparison with the CMRs of lower redshift clusters 
shows that the red galaxy overdensity primarily lies blue- 
ward of the no-evolution predictions. That the red galax- 
ies in 1138 are also redder than the 

^•forrn ^ 3 case sug- 
gests both that there are galaxies with significant dust 
content, an assertion supported by the SED fits, and also 
that they were perhaps formed at higher redshift than the 
RDCS1252 galaxies. Of course, without a classical, low- 
scatter red-sequence to use as a baseline there remains 
considerable uncertainty in the age of the population as 
a whole. The results of Steidel et al. (2005) suggest that 
protocluster galaxies are older than their "field" counter- 
parts at z ~ 2.3 and that these ages and stellar masses 
were broadly consistent with evolution to lower redshift 
cluster galaxies. However, their protocluster members 
were all UV-selected and star-forming. With future spec- 
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TABLE 1 

Properties of Red Galaxies in the MRC 1138-262 Field: Hiqq < 24.5 A 1.1 < (./no — Hl6o) < 21 
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(AB) 
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0.98 


1.00 


3.4 ± 
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1.79 ± 
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0.24 
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0.01 


1.99 ± 
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22.50 ± 
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28.88 


-27.16 
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1.52 


314 
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1.00 


1.00 
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0.7 


0.15 
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0.01 


1.28 ± 





05 


23.03 ± 


01 


1.92 


1.38 


1.89 


1.29 
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2.10 


1.00 


1.00 
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0.3 


0.40 


± 


0.03 


3.34 ± 





23 


21.77 ± 0- 


01 


1.80 


1.89 


2.22 


1.61 
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2.23 
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1.00 
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0.12 


± 


0.00 


1.02 ± 





04 
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01 


1.30 


1.13 


1.67 


1.36 
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2.25 


1.00 


1.00 


3.5 ± 


0.3 


0.56 


± 


0.04 


4.72 ± 
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22.75 ± 0. 


02 


0.63 


1.28 


2.19 


1.64 
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2.42 


1.00 


1.00 


5.0 ± 
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0.38 
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0.02 


3.20 ± 
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22.00 ± 0- 


01 


-28.02 


1.45 
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2.44 


1.00 


1.00 


2.7 ± 
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0.12 
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0.02 


1.01 ± 
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22.64 ± 0- 
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1.04 


1.37 
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1: quiescent Total magnitudes based on the SExtractor MAG_AUTO pa 



troscopy of our red galaxy sample it will be possible to 
see if these differences persist when looking at a more 
varied galaxy sample. 

For three samples of galaxies drawn from the full H\qo- 
band selected dataset we have fit a color-magnitude re- 
lation and estimated the intrinsic scatter. The CMR at 
z = 2.16 is not as well-defined as at z ~ 1 or 0. For 
sample C, made up of 8 galaxies, the color, best-fit spec- 
tral template, morphology and photo- z all point towards 
them being (proto-)clliptical galaxies within the proto- 
cluster. For this small sample, the estimated intrinsic 
scatter is rather low and may suggest that these galax- 
ies represent the forming red-sequence in this protoclus- 
ter. The slope of this relation is extremely steep com- 
pared to lower redshift clusters. The slope of the CMR 
is generally assumed to be a manifestation of the mass- 
metallicity relation and would therefore flatten at higher 
redshift. The major caveat regarding the steep slope of 
Sample C is that none of these galaxies are spectroscop- 
ically confirmed protocluster members. Therefore, this 
"relation" may just be a random, although somewhat 



unlikely, coincidence rather than a nascent CMR. How- 
ever, further deep NIR imaging coverage of this field is 
required to identify additional members of this proto- 
elliptical galaxy class. 
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